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Abstract. We consider here the constraints in SUGRA models on the SUSY parameter space 
due to current experimental bounds on the light Higgs mass m/j, the b — > sj decay, the amount 
of neutralino cold dark matter ilh 2 , and the muon magnetic moment. Models with universal 
soft breaking (mSUGRA) and non-universal gaugino or Higgs masses at Mq are examined. 
For mSUGRA, the mu, b — > sj and ilh 2 constraints imply a lower bound on the gaugino mass 

of rax/2 ^300GeV implying the gluino and squarks have mass ~700GeV, and the neutralino 

~120GeV. The current status of the Brookhaven muon g - 2 experiment is reviewed, and if the 
Standard Model (SM) contribution evaluated using the e + + e~ data is correct, a 2a bound on 
the deviation of experiment from the SM produces an upper bound on m 1( i 2 that eliminates the 
"focus point" regions of parameter space. Dark matter (DM) detection cross sections range 
from 5 x 10~ 8 pb to 5 x 10~ 10 pb which would be accessible to future planned detectors. The 
SUSY decay B s — > /i + + \i~ is seen to be accessible to the Tevatron Run 2B with 15 fb _1 

luminosity for tanbeta ~ 30. The most favorable signals of SUSY for linear colliders are stau 
pair production and neutralino pair production, though it will require an 800GeV machine to 
cover the full parameter space. Non-universal models can modify some of the above results. 
Thus a non-universal (heavier) gluino mass at Mq can significantly reduce the lower bound 
constraints of b — > sj and rrih giving rise to a lighter SUSY spectrum. A heavier up Higgs 
mass can open an additional region with allowed relic density arising from annihilation via the 
s-channel Z diagram with an O(10) larger DM detector cross section. 



1. Introduction 

We are now approaching the time when experiments will tell us what the new physics is that 
will replace the Standard Model. Neutrino experiments have already shown the breakdown of 
the Standard Model with the evidence of neutrino masses. However, this does not shed light on 
how the new dynamics resolves the gauge hierarchy problem, and there are many suggestions 
for this: SUSY and SUGRA models, extra dimensions (large and small), string/M-theory 
models, etc. 

SUSY models all contain a large number of parameters to be determined finally by 
experiment. The more phenomena the model applies to, the more one will be able to limit 
the parameter space. We consider here SUSY models with R-parity invariance, since these 
have the unexpected prediction of the existence of dark matter [[l], §]. The requirement that 
the theory predict the experimentally measured amount of dark matter already puts strong 
constraints on the parameter space. The MSSM with over 100 free parameters is not very 
restrictive. At the other extreme is mSUGRA which is consistent with grand unification (up 
to now the only direct experimental evidence in favor of SUSY!) and has only four new 
parameters and one sign. It's natural therefore to start with mSUGRA and see how the current 
experiments restrict its parameter space, and hence what predictions can be made for future 
dark matter experiments and for accelerator experiments at the LHC and NLC. We will then 
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perturb the mSUGRA framework with non-universal soft breakings that are consistent with 
current data, and thus get some idea of how robust the mSUGRA predictions are. 

The current experiments that most strongly restrict the SUSY parameter space are the 
following: 

(1) The amount of cold dark matter (CDM). Global fits to the CMB and other data yield 
^cdm^ 2 = 0.139 ± 0.026 0], and we take a 2.5a range around the central value: 

0.07 < ttcDMh 2 < 0.21 (1) 

The MAP data, due out in the near future, will be able to significantly narrow this range. 

(2) Higgs mass. The LEP lower bound of nih > 114.1 |Q] is a significant constraint for 
lower tan/3 (i. e. tan/5 ~ 30), but generally will become very significant if the bound were 
to rise by just a few GeV. Unfortunately, theoretical calculations still have an error of about (2 
- 3) GeV, and so we will (conservatively) interpret this experimental bound to mean that that 
the theory calculation of m h should obey m h > 111 GeV. mSUGRA predicts m h ~ 130 GeV, 
which could make the light Higgs within the reach of the Tevatron Run2B. 

(3) b — ► S7 decay. The CLEO data[|]] has both systematic and theoretical error, and so 
we take a relatively broad range around the central value: 

1.8 x 1(T 4 < B(B -> X s7 ) < 4.5 x 10~ 4 (2) 

The b — > S7 rate is significant for large tan f3 (tan (3 ~ 30) where it produces a lower bound 
on the gaugino mass m\/2- If the lower bound on the branching ratio were raised, it would 
increase this lower bound on m\/2- 

(4) Muon magnetic moment anomaly, a^. Since this conference, the Brookhaven E821 
experiment has released new results that have halved the statistical error with the central 
value essentially unchanged [§]. Also, new data from CMD-2 and BES has allowed a better 
evaluation of the SM contribution and two new theoretical evaluations have appeared [|7|, §]. 
Using the e + — e~ data to evaluate the hadronic contribution, both analyses lead now to a 3a 
deviation, e.g. [[7]]: 

Aa M = 33.9(11.2) x 1(T 10 (3) 

If, on the other hand, the tau decay data (with CVC breaking corrections included) is used, 
one finds only a 1.6a discrepancy [[7J, and further, the two results are statistically inconsistent. 
In the following, we will assume that the e + — e~ data is correct, and take a 2a band around 
the central value of Eq.(3) to see what this implies. However, it may turn out that the tau data 
evaluation is correct, under which circumstances, the current a /t data will give no significant 
constraint. 

Future experiments at accelerators will eventually determine if SUSY is correct. Thus 
Run2B at the Tevatron should be sensitive to B s — > fi + + fi~ with 15fb~Vdetector data [Q, 
and perhaps distinguish different SUGRA mediation models [JTOt] . The LHC can, of course 
measure SUSY masses and couplings, e.g. for the gluino up to 2.5 TeV (i.e. ~ 1 TeV), 
while the NLC might observe stau and neutralino production processes. 

2. mSUGRA Model 

The mSUGRA model depends upon four parameters and one sign. It is convenient to chose 
these as follows:(i) m , the scalar soft breaking mass at the GUT scale Mq- (2) m^, the 
gaugino mass at M G . [Note that m^o = 0.4 m 1 / 2 (where x\ ls m e lightest neutralino), 
m~± = 0.8mi/2 (where Xi 1S me lightest chargino), and m g = 2.5mi/ 2 (where g is the 
gluino).] (3) A the cubic soft breaking mass at M G . (4) tan (3 =< H 2 > / < Hx > at the 



electroweak scale, where if 2 ,i gives mass to the (up,down) quarks. In addition there is the 
sign of //, the Higgs mixing parameter in the superpotential term W = fxHiH 2 . We examine 
the parameter space over the following range: m > 0, my 2 < ITeV (m§ < 2.5 TeV), 
2 < tan/3 < 55, and \ A Q \ < 4my 2 . 

For heavy nuclei, the spin independent neutralino-nucleus cross section dominates for 
terrestrial dark matter detectors, which allows one to extract the Xx - proton cross section 
CyO-p- (Here xl ls me dark matter candidate.) The neutralino scattering by quarks in the 
nuclear target proceeds through s-channel squark states, and t-channel Higgs boson states 
(h, H). To calculate the relic density of neutralinos left over after the Big Bang, one needs 
the neutralino annihilation amplitudes. The annihilation proceeds through the s-channel Z 
and Higgs channels (h, H, A), and t-channel sfermion channels. However, if a second 
particle is nearly degenerate with the Xx, one must include it in the early universe annihilation 
processes, which leads to the co-annihilation phenomena. In SUGRA models this accidental 
near degeneracy occurs naturally for the light stau, f 1 . Co-annihilation then begins at my 2 = 
(350-400)GeV and the scalar mass m must be raised as mx/ 2 increases to keep the fx heavier 
than the xl- 

The neutralino -proton cross section has the following general behavior: cr^o_ p increases 
with increasing tan (3, and decreases with increasing mx/2, m o- The maximum value of cr^o-p 
then generally occurs at large tan (3 and small mx/2, m o- 

One starts the analysis at M G and uses the renormalization group equations (RGE) to 
obtain predictions at the electroweak scale. In carrying out these calculation, it is necessary 
to include a number of corrections and we list some of these here: (1) We use two loop gauge 
and one loop Yukawa RGE in running from M G to the electroweak scale M EW , and three loop 
QCD RGE below Mew for light quark contributions. (2) Two loop and pole mass corrections 
are included in the calculation of m h . (3) One loop correction to m b and m T are included [pTj]. 
(5) All stau-neutralino co-annihilation channels are included in the relic density calculation 
[ |T2"| , [TBI , (Chargino-neutralino co-annihilation does not occur for mx/2 < 1 TeV.) Large 
tan/3 NLO SUSY corrections to b — > S7 are included [|T5|, [T5p. We do not include Yukawa 
unifications or proton decay constraints as these depend sensitively on post-GUT physics, 
about which little is known. 

Figs. 1 and 2 illustrate the allowed regions in the m — mx/2 plane for tan f3 — 10 and 
tan/3 = 50 for A = 0, /i > 0. We note for tan/3 = 10, that the lower bound on my 2 is 
set by the lower bound on the Higgs mass. The narrowness of the dark matter allowed band 
does not imply any fine tuning. Thus the lower edge is determined from the co-annihilation 
effect so that the fx, which is nearly degenerate with the x\ does not cause too much early 
universe annihilation (violating the lower bound of Eq(l)), while above the upper bound of 
the band, too little early universe annihilation occurs (violating the upper bound of Eq. (1)). 
Thus the MAP (and Planck) satellite data will narrow this band further, with more accurate 
determinations of ^cdm- At tan/3 = 50, the allowed dark matter band expands at low mx/ 2 
due to the fact that the A Higgs becomes light, allowing additional early universe annihilation, 
and this is compensated by raising m . Note here it is the b — > 37 branching ratio that give 
the lower cut off in mx/ 2 . The vertical lines through the dark matter allowed bands represent 
the values of cr^,o_ p . We see that the cross sections range from 5 x 10~ 8 pb to 10 -9 pb. The 
fi < possibility (not shown) allows for regions of much smaller dark matter cross sections 
QT7| , |T2p. However, if the SUSY contribution to a M is indeed positive (which is the case for the 
e + — e~ evaluation of the SM contribution), then the ji < possibility is eliminated [fi~8L |19|]. 
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Figure 1. Allowed regions in mO — mu2 parameter space for tan/3 = 10, A$ — 0, \i > 0. 
The left vertical line corresponds to c^o_ p = 5 x 1CT 9 pb and the right to 1CT 9 pb. 



3. The Muon Magnetic Moment Anomaly 

As discussed in the introduction, it is still unclear if the current measurements of the muon 
magnetic moment indicate a deviation with the Standard Model prediction. In this section 
we will assume that the analysis based on the e + — e~ data, which yields the 3a deviation of 
Eq.(3), is correct, in order to see what restriction on the SUSY parameter space this possibility 
allows. In supersymmetry, one expects a contribution to the muon magnetic moment of 
characteristic size 10 x 1CT 10 arising from chargino-sneutrino and neutralino-smuon loops 
(with the photon attached to any charged line). If the final result requires a^ USY to be much 
smaller than 10 x 10~ 10 , it would mean that the SUSY squark and gluino masses are in the 
TeV domain, while a value of a^ USY ~ 40 x 10~ 10 would exclude mSUGRA [0. Thus the 
allowable range for a^ USY is already somewhat restricted. Fig. 3 exhibits the effect of the muon 
magnetic moment anomaly for the case of tan/3 = 40, A = 0, /x > 0. The shaded region in 
the upper right would be the part of the parameter space excluded at the 2a level (assuming 
the e + — evaluation of the SM). Combined with the dark matter constraint, one sees that the 
constraint give rise to an upper bound mi/2 ~ 800GeV or m^o < 320 GeV, m~± < 640 
GeV and m- g < 2 TeV. The b — > 57 constraint gives a lower bound of > 350GeV, 
and thus the parameter space is constrained at both ends. Note also that it would not take a 
significantly larger value then the central value of 33 x 10 -10 to exclude the full parameter 
space for this value of tan (3. Thus a resolution of the Standard Model contribution to is of 
great importance for SUSY models. 
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Figure 2. Same as Fig. 1 with tan j3 = 50, and cr^o_ p = 5 x 10~ 8 (left) and 2 x 1CT 9 (right). 



4. The B s -> //+ + \r Decay 



The 5 S — > fi + + fi~ decay pT| , ^, offers another window for investigating the 

mSUGRA parameter space. The Standard Model predicts a branching ratio that is quite small 
i. e. B[B S — > fi + + = (3.1 ± 1.4) x 1CT 9 . However, the SUSY contribution can become 
quite large for large tan (3. An example of one of the important graphs for this process is 
shown in Fig. 4, where the amplitude grows as (tan /3) 3 , and hence the rate grows as (tan (3) 6 . 
This process has become interesting because it appears possible to observe it at Run2B of the 
Tevatron [j^]. A set of cuts eliminating the background (e.g. gluon splitting of g — » bb and 
non-b background) exits and leads to a sensitivity for CDF for a luminosity of 15fb — 1 of 
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(and if a similar sensitivity can be obtained by DO the combined sensitivity would be 
0.65 x 10~ 8 ). Fig. 5 shows the branching ratio for different tan [3 for the case m = 300GeV, 
A = 0, fi > 0, and Fig. 6 shows the expected limit on the branching ratio observable by CDF 
as a function of luminosity. We see that with 15 fb _1 , CDF should be sensitive to this decay 
for tan (3 ~ 30, and in fact too large a branching ratio could exclude mSUGRA with 2 fb _1 . 



5. NLC Reach 

We consider here two possible accelerators at energies of 500GeV and 800GeV to examine 
what part of the SUSY spectrum would be available to linear colliders. We've seen that 
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Figure 3. Allowed region in the mo — mi/% plane for tan (3 — 40, A$ — 0, fi > 0. The shaded 
upper region is excluded at the 2er level for a M obeying Eq.(3). 




Figure 4. A leading diagram for the decay B s — > + /i . The vertices with the heavy dot 
are proportional to tan (5. 



the and 6 — > S7 bounds already means that ~ (350 — 400) GeV, and so for 

mSUGRA, gluinos and squarks would generally be beyond the reach of such machines. The 
most favorable SUSY signals are then 

e + + e" - x°2 + Xi - + I' + Xi) + Xi (5) 

e + + e" -> f+ + ff -> (r + Xi) + (r + X?) 
where Z 1 * 1 is any charged lepton. Since in mSUGRA m^o ~ 2m^o one has for the kinematic 
mass reach l/2m^o ~ m^o ~ 165(265) GeV for ^/~s = 500(800) GeV, and m fl ~ 250(400) 
GeV for y/s = 500(800) GeV In general the stau is the lightest slpeton and so will be the 
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Figure 5. The branching ratio of B s 
tan/3[9]. 



fi as a function of mx/2 for different values of 



major signal, as the selectron and smuon will not be accessible for most of the parameter 
space. There are a number of backgrounds however, and appropriate cuts have to be arranged 
to suppress them. This is currently being studied [bdutta4]. 

6. The mSUGRA Parameter Space 

We now summarize the effects on the mSUGRA parameter space of the constraints from mu, 
b — > sj, dark matter, and (assuming Eq. (3) is valid). We will also see what effect might 
be expected from measurement of the B s — > fi/i decay at the Tevatron, and examine the 
NLC reach for SUSY. Fig. 7 exhibits these effects for tan (3 =10, A = 0, \i > 0. We see 
that for low tan (3, the constraint eliminates most of the parameter space, and the NLC 
with y/s = 500 GeV would be able to scan the full remaining space with either the stau or 
neutralino signal. The <J^o_ p cross sections are small, but of the size that future experiments 
hope to achieve (e.g. GENIUS, Cryoarray, ZEPLIN IV, CUORE).Of the two NLC signals, 
the ~~ Xi i s sensitive to large m , while the f\ — f\ is sensitive to relatively large m\/2- Due 
to the dark matter constraint, the latter appears to be the more important signal. 

Fig. 8 shows similar information for tan (3 = 40. For larger tan (3, the possibility 
of observing the B s — > fxfi signal at the Tevatron becomes significant, and assuming the 
combined CDF and DO data were available, it could cover the full parameter space for an 
anomaly > 10 x 10~ 10 . The NLC (500GeV) tau-tau signal would now cover only about 
one half the allowed parameter space, while for 800GeV the full parameter space could be 
examined. 




For very large tanbeta, a new feature of a "bulge" occurs in the dark matter channel at 
low m 1/2 in the allowed dark matter channel. This is due to the fact the heavy Higgs A and H 
become light, and allow a more rapid annihilation in the early universe through the A and H 
s-channel poles. Figs. 9 and 10 for tan (3 = 50 and 55 exhibit this feature. (We note in general 
that the results for tan (3 > 45 become very sensitive to the precise values of m t and m^. We 
are here using the central values of m t (pole) = 175 GeV and m b {m b ) = 4.25GeV, but the 
figures would change significantly if one were to go one or two a away from the mean.) We 
see from the figures that the fi — fx LC signal would be able to cover the full parameter space 
of tan (3 = 50 for an 800 GeV collider, but not quite cover it all for tan (3 = 55. Of course, 
for these large tanbeta, the B s — > /i/i signal is greatly enhanced, as is the dark matter detector 
cross sections. 



7. Non-universal Models 

Our previous discussion has been within the framework of mSUGRA with universal soft 
breaking occurring. It is interesting to investigate what aspects of these results would survive 
in the presence of some non-universal soft breaking. One peculiar feature of mSUGRA is 
the narrow dark matter allowed bands due in part to the "accidental" near degeneracy of the 
fi and the Xi- This phenomena would be maintained even if the gaugino masses were non- 
universal at M G , since m\ — m 2 ~ depends mainly on rh\. Thus this type of co-annihilation 

effect is generic. A second feature of mSUGRA, that the combined effects of the m h and 
b — > 57 experimental bounds put a lower limit m\/2 > (300 — 400) GeV for all tanbeta, is 
however sensitive to the mSUGRA assumption of universal gaugino masses at M G . Thus if 
one assumes gluino breaking of this universality at M G , 

rhi = fa 2 = m 3 (l + 5 3 ), (6) 

then a positive 5 3 raises the stop mass and effects strikingly both the b — > S7 branching ratio 
and the value of m h . This is illustrated in Fig. 1 1 for the case of <5 3 = 1, tan (3 = 50, A = 0, 
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Figure 7. Allowed region in the too — mi /% plane for tan (3 — 10, A$ = 0, /i > 0. The shaded 
upper region is the 2a bound. The vertical dash-dot lines are the kinematical bound at a 
linear collider for the Xi ~ x\ signal, the left line for a 500 GeV machine and the right for 800 
GeV. (Parameter space to left of these lines would be observable if the kinematical reach could 
be achieved.) The curved lines are similarly for the T\ — f\ signal at a linear collider. The solid 
vertical lines through the dark matter allowed band give the values of <J^a_ p i.e. 5 x 10 _9 pb 

for the left line and 1 x 10~ 9 for the right line. 



H > 0. One sees that both the b — > s^y and m h constraint is moved strongly towards lower 
mi/2, and now the lower bound on mi/ 2 is only mi/ 2 ~ 190 GeV. The Higgs lines are also 
moved and lie approximately 150 GeV lower in m^. The a M bound, however, becomes 
somewhat more constraining. 

A second non-universality that might naturally arise is with the scalar soft breaking 
masses. Thus while flavor changing neutral current constraint require squark masses to be 
nearly degenerate at the GUT scale, there is little theoretical reason to require that the Higgs 
masses be degenerate. We consider here the simple model where at Mq 



where m is the universal squark and slepton masses. While this introduces two new 
parameters into the model, one may qualitatively understand the effects of 5i >2 . In SUGRA 
models, the fi parameter governs much of the physics, and if /j 2 decreases (increases), then 
the Higgsino content of the neutralino xi increases (decreases). The RGE with radiative 
breaking shows that /i 2 is indeed sensitive to <5 2 (but only slightly sensitive to 5±). For # 2 
positive, /j 2 decreases, and this produces two effects: (1) it increases the dark matter detection 
cross section cr%o_ p since this quantity depends on the interference between the Higgsino and 
gaugino parts of the neutralino. (2) It increases the x^ — x^ — Z coupling, allowing more early 
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Figure 8. Same as Fig. 7 with dashed lines giving B s — > /i/i Tevatron reaches. The left solid 
dark line corresponds to o^-p = 3 x 1CP 8 and the extreme right to 1 x 10~ 9 . 
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Figure 9. Same as Fig. 8 for tan/3 = 50, with the left solid line corresponding to 
cr^o_ p = 5 x 10~ 8 pb and the right to 2 x lCT 8 pb. 
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Figure 10. Same as Fig. 8 for tan j3 = 55. 



universe annihilation of the neutralinos through the Z s-channel. This then opens another 
channel of allowed relic dark matter in the m — mi/ 2 plane at relatively low mi/ 2 but high 
m (so that too much annihilation doesn't occur). These effects are seen in Figs. 12 and 13. 
Fig. 12 shows the new Z-channel allowed dark matter region at high m , along with the usual 
stau co-annihilation narrow band at low m for tan (3 = 40. Fig. 13 shows the corresponding 
values of cr^?~p' me u PP er dashed band corresponding to the Z-channel of Fig. 13. One sees 
that the detector cross section for this possibility can be a factor of ten or more larger than the 
usual stau band. 



8. Conclusions and Summary 

We have examined here for SUGRA models of neutralino dark matter, the constraints on the 
SUSY parameter space that arise from current experiments, and have considered models both 
with universal soft breaking (mSUGRA) and with non-universal soft breaking. 

For mSUGRA, co-annihilation effects and the current bounds on the relic density 
generally restrict the allowed parameters to be in a narrow band in the m — mi/ 2 plane 
with relatively low m (except in the low m 1 /2 region where for very large tan/5 a "bulge" 

can exist). If the SUSY contribution to the muon magnetic moment anomaly is ~ 10 x 10~ 10 
(which is the current 2a bound for the e + — e~ data analysis of the SM contribution^, |8p), 
then the "focus point" region [ p5| ] of very large m is also eliminated. The current bounds 
on rrih and b — > further produce a lower bound on mi/ 2 of mi/ 2 ~ (300- 350) GeV over 
the full tanbeta domain. This implies that m^o ~120 GeV, M-± ~ 240 GeV and m~ g ~ 750 




GeV. These mass bounds are larger than what exists now from accelerators. In general we 
find that dark matter direct detection neutralino -proton cross sections range from ~5x 10~ 8 
to ~ 5 x 10 10 , which are mostly within the range of sensitivity of future planned dark matter 
detectors. The B s — > decay will be accessible to the Tevatron over much of the parameter 

space for tanbeta ~ 30 with 15 fb _1 of luminosity. The accessibility of the f\ — f\ and Xi — X2 
SUSY signals for the NLC was examined. 

To see how robust the above results are, we examined two types of non-universal soft 
breaking models: those with non-universal gaugino masses at the GUT scale, and those with 
non-universal Higgs masses. In general for these types of models, the co-annihilation effects 
still remain. However, an increased non-universal gluino mass can significantly reduce the 
lower bound constraints on mi/ 2 arising from the b — > S7 and rrih constraints, making these 
constraints less significant. ( The a M constraint is simultaneously somewhat strengthened). 
Non-universal Higgs soft breaking masses can also produce some striking effects. Thus an 
increase of the H 2 mass at the GUT scale increases the Higgsino part of the neutralino, which 
then increases the x\ ~ Xi ~ Z coupling and opens a new Z s-channel annihilation region at 
low mi/2 and high m with acceptable amount of relic dark matter. In this region the Xi ~ V 
detection cross section is a factor of 10 or more larger than that from the fi co-annihilation 
region, making it accessible to the next round of dark matter detector experiments. However, 
SUSY signals in such regions would be difficult to detect for a 500 GeV NLC. 
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Figure 12. Allowed region in the mo — mi/2 plane for 82 = 1, tan/3 = 40, Ao = mxii, 
/i > 0. The lower narrow band is the ususal stau co-annihilation band, and the upper band is 
due to the Z-channel annihilation. 




Figure 13. Dark matter detection cross section for the two allowed regions in the mo — m ly i 2 
of Fig. 12. The upper dashed band is for the Z-cannel annihilation, and the lower line is for 
the fi — X? co-annihilation channel. [12] 
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